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Abstract. It is well known that for many semilinear parabolic equations there 
is a global attractor which has a cell complex structure with finite dimensional 
cells. Additionally, many semilinear parabolic equations have equilibria with 
finite dimensional unstable manifolds. In this article, these results are unified 
to show that for a specific parabolic equation on an unbounded domain, the 
space of heteroclinic orbits has a cell complex structure with finite dimensional 
cells. The result depends crucially on the choice of spatial dimension and the 
degree of the nonlinearity in the parabolic equation, and thereby requires some 
delicate treatment. 



1. Introduction 
In this article, the space of heterochnic orbits of 

du{t,x) d'^u(t,x) y, . ,, , 

is shown to have the structure of a cell complex with finite-dimensional cells, where 
u G Ci(K,C°'"(R)), (j) € L'^n C°'"(R), and |0| ^ as |a:| ^ oo. This result is a 
generalization of a well-known result that the unstable manifolds of ([T|) are finite 
dimensional. Indeed Theorem 5.2.1 in [TD] can easily be made to apply with the Ba- 
nach spaces we shall choose. Theorem 5.2.1 in [Tn| shows the existence of a smooth 
finite dimensional unstable manifold locally at an equilibrium. One can then use 
the iterated time-1 map of the flow for ([1]) to extend this local manifold to a max- 
imal unstable manifold. There are also finite Hausdorff dimensional attractors for 
the forward Cauchy problem on bounded domains [16] . However, we shall exhibit 
a more global approach to the finite dimensionality of the unstable manifolds. This 
approach allows us to examine the finite dimensionality of the space of heteroclinic 
orbits connecting a pair of equilibria, which is a new result in the spirit of [5] . The 
techniques used here depend rather delicately on both the degree of the nonlinear- 
ity (which is quadratic) and the spatial dimension (which is 1). Both of these are 
important in the standard methodology as well, as the portion of the spectrum of 
the linearization in the right half-plane needs to be bounded away from zero. In 
the case of ([1]), the spectrum in the right-half plane is discrete and consists of a 
finite number of points. 

Of an immediate and important concern is that there may not be any solutions 
to ll]) which are defined in C^(M, C"'"(M)). More particularly, are there solutions to 
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([T]) which are defined for aU time? This question can be answered in the affirmative 
[15j , so this article makes the assumption that the space of heterochnes is nonempty. 

2. Applications 

Equation ([T]) is a very simple model of combustion. If is a positive constant, 
then the equation supports traveling waves. Such traveling waves can model the 
propagation of a flame through a fuel source pTl . 

In addition to a model of combustion, ([T]) can also be a simple model of the 
population of a single species, with a spatially- varying carrying capacity, 4>. Indeed, 
one easily finds that under certain conditions the behavior of solutions to ([1]) is 
reminiscent of the growth and (admittedly tenuous) control of invasive species [2]. 
It is the control of invasive species that is of most interest, and it is also what the 
structure of the attaching maps of the cell complex reveals. In the example given in 
Section [6l there is one more stable equilibrium, and several other less stable ones. 
The more stable equilibrium can be thought of as the situation where an invasive 
species dominates. The task, then, is to try to perturb the system so that it no 
longer is attracted to that equilibrium. An optimal control approach is to perturb 
the system so that it barely crosses the boundary of the stable manifold of the the 
undesired equilibrium, and thereby the invasive species is eventually brought under 
control with minimal disturbance to the rest of the environment. 

3. Prior work 

Equations of the form ([1]) have been of interest to researchers for quite some 
time. Existence and uniqueness of solutions on short time intervals (on strips 
(— to,io) X I^) can been shown using semigroup methods and are entirely standard 
[21]. However, there are obstructions to the existence of eternal solutions. Aside 
from the typical loss of regularity due to solving the backwards heat equation, there 
is also a blow-up phenomenon which can spoil existence in the forward-time solution 
to ([T]) . Blow-up phenonmena in the forward time Cauchy problem (where one does 
not consider t < 0) have been studied by a number of authors [9] [6] [22] [13] [3] 
[26] [27] . More recently, Zhang et al. ([25] [20] [23]) studied global existence for the 
forward Cauchy problem for 

^ ^ Au + u^ - V(x)u 
at 

for positive u, V . Du and Ma studied a related problem in [S] under more restricted 
conditions on the coefficients but they obtained stronger existence results. In fact, 
they found that all of the solutions which were defined for all t > tended to 
equilibrium solutions. 

The boundary value problem that results from taking a; G C M" for some 
bounded (instead of a; £ M") has also been discussed extensively in the literature 
[TU] [TT] [3]. For the boundary value problem, all bounded forward Cauchy problem 
solutions tend to limits as \t\ oo, and these limits are equilibrium solutions. 

Almost all of the literature (including this article) describing eternal solutions 
to ([I} is restricted to discussing heterochnes. For unbounded domains and certain 
symmetries of 0, one can find traveling waves. Since the propagation of waves in 
nonlinear models is of great interest in applications, there is much written on the 
subject. The general idea is that one makes a change of variables (i, x) i— > ^ — x — ct 
which reduces ([1]) to an ordinary differential equation. This ordinary differential 
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equation describes the profile of a traveling wave. Powerful topologically-motivated 
techniques, such as the Leray-Schauder degree, can be used to prove existence of 
wave solutions to ([T|). Asymptotic methods can be used to determine the wave 
speed c, which is often of interest in applications. See [21] for a very thorough 
introduction to the subject of traveling waves in ([T]). 

4. The linearization and its kernel 

We begin by considering an equilibrium solution / to ([T]). As discussed in [17], 
this solution has asymptotic behavior which places it in fl n L°°{M.). We are 
particularly interested in solutions which lie in the a-limit set of /, those solutions 
which are defined for alH < and tend to /. Center attention on this equilibrium 
by applying the change of variables u{t, x) i-^ u{t, x) — f{x) to obtain 

' §-^u{t, x) = ^u{t, x) ~ 2f{x)u{t, x) - u^{t, x) 
^ u{0,x) = h{x) e C2(M) 
limt^_oo u(t, x) = f{x) 
< 0,a; e M. 

Thus we have a final value problem for our nonlinear equation. All solutions to 
(12) will tend to zero as t — > — oo uniformly by Lemma 6 of |TS]. Of course, ^ is 
ill-posed. We show that there is only a finite dimensional manifold of choices of h 
for which a solution exists. 

4.1. Backward time decay. The decay of solutions to zero is a crucial part of the 
analysis, as it provides the ability to perform Laplace transforms. In the forward 
time direction, one obtains upper bounds for solutions by way of maximum princi- 
ples, and lower bounds for the upper bounds by way of Harnack estimates. In the 
backward time direction, these tools reverse roles. Harnack estimates provide upper 
bounds, while the maximum principle provides lower bounds for the upper bound. 
In the proof of Lemma 6 of [18], the latter was used to some advantage. In this 
section, we briefly apply a standard Harnack estimate to obtain an exponentially 
decaying upper bound. 

Harnack estimates for a very general class of parabolic equations are discussed 
in [1^ and [1]. In those articles, the authors examine positive solutions to 

div A{x, t, u, "S/u) — — — B{x, t, u, Vu), 

where a; G M", and A : M2"+2 ^ M« and B : M2"+2 ^ k satisfy 
\A{x,t,u,p)\ < a\p\ + c\u\ + e 
\B{x,t,u,p)\ < b\p\+d\u\+f 

p-A{x,t,u,p) > 

a 



for some a > and b, ...g are measurable functions. For a solution u defined on 
a rectangle R, the authors define a pair of congruent, disjoint closed rectangles 
C B with being a backward time translation of The main result 
is the Harnack inequality 

(3) maxu < 7 ( minu + L \ , 
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where 7 > depends only on geometry and a (but not b, ...g) and L is a linear 
combination oi e, f, g whose coefficients depend on geometry. 

In the case of (H]), or indeed of the analogous equation with higher degree terms, 
we have that ([3]) will apply with L — 0. Notice that the conditions on A, B are 
satisfied because any solution to ([2]) is automatically a finite energy solution, and 
therefore is bounded and has bounded first derivatives. The only difficulty is that 
([3]) applies for positive solutions, while ([2|) may have solutions with negative por- 
tions. However, one can pose the problem for the (weak) solution of 

sgn (u) (Au - - 2/u) 

A|u| -2/|u| 
AH-|up-2|/|H 



dt 



> 



for which we only get positive solutions. By iterating ([3]) we have that solutions to 
([2]) decay exponentially as t — > — cx). 

4.2. Topological considerations. 

Definition 1. Let Ya{X) be the subspace of C^{X,C°^°'{R)) which consists of 
functions which decay exponentially to zero like e"*, where < a < 1. We define 
the weighted norm 

My. = ||e-'^1w(i)llco.»(M)||ci 

and the space 

YaiX) = {u = u{t,x) e C'iX,C"-"{R))\\\u\\Y. < 00} . 

In a similar way, we can define the weighted Banach space Za{X) as a subspace of 
C°(X,C°'"(]R)). It is quite important that Ya and Za are Banach algebras under 
pointwise multiplication. 

In light of the previous section, solutions to ^ are zeros of the densely defined 
nonlinear operator N : ya((— oo,0]) — > Za{{— 00, 0]) given by 

(4) + 

About the zero function, the linearization of N is the densely defined linear map 
L : Fa ((-00, 0]) -> Z,((-oo, 0]) given by 

d d'^ d 

where we define H = ^3 — 2/. Also note that L is the Frechet derivative of N, 
which follows from the fact that Ya and Za are Banach algebras. 

Remark 2. We are using C°'"(R) instead of C°(M) to ensure that N and L be 
densely defined. We could use space of continous functions which decay to zero, or 
the space of uniformly continous functions equally well. 



Convention 3. We shall conventionally take a > to be smaller than the smallest 
eigenvalue of H . 
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We show two things: that the kernel of L is finite dimensional, and that L 
is surjective. These two facts enable us to use the implicit function theorem to 
conclude that the space of solutions comprising the a-limit set of an equilibrium is 
a finite dimensional submanifold of ya((— oo, 0]). 

4.3. Dimension of the kernel. 

Lemma 4. // / is an equilibrium solution, then the operator L : Ya{{~oo,0\) — > 
00,0]) in ([5|) has a finite dimensional kernel. 

Proof. Notice that the operator L is separable, so we try the usual separation 
h{t,x) = T{t)X{x). Substituting into jS]) gives 

f d \ 



T ~ X 



= A 



for some A G C. The separated equation for T yields T = CajC^*. Since we are 
looking for the kernel of L in Yq C L°°(IR^), we must conclude that A must have 

nonnegative real part. On the other hand, the spectrum of H = (^-^ — 2/^ is 

strictly real, so A > 0. Indeed, there are finitely many positive possibilities for A each 
with finite-dimensional eigenspace. This is a standard fact about the Schrodinger 
operator H since / is an equilibrium. Thus L has a finite dimensional kernel. □ 

4.4. Surjectivity of the Unearization. In order to show the surjectivity of L, 
we will construct a map T : oo,0]) yQ((— oo,0]) for which LoT = idz^. 

That is, we construct a right-inverse to L, noting of course that L is typically not 
injective. We shall derive a formula for T using the Laplace transform u i— > ?7 

v{s,x)— / e''*v{t,x)dt, 



(6) 



where ^{s) > —a and v g oo,0]). 

Since Lemma 2] essentially solves ^ , we will be solving the inhomogeneous 
problem with zero final condition 

I T - + 2/ ^) - ^w{t, x) e ZM-oo, 0]) 

\w(0,a;) = 

for t < 0. The Laplace transform of this problem is 

sv{s,x)-\ ^ 2f{x)v{s,x) = w{s,x) 

{H + s)v{s,x) = w{s,x). 

Choose a vertical contour C with > K(s) > —a, so that the Laplace transforms 
are well-defined, and that the contour remains entirely in the resolvent set of —H. 
Then we can invert to obtain 

v{s, x) = {H + s)^^w{s, x). 
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Figure 1 . Definition of the contour C 

Using the inversion formula for the Laplace transform yields 

v{t, x) 



— I e-'\H + s)-^w{s,x)ds 
2vri Jc 

— I e"\H + s)-^ I e'^w{T,x)dTds 
Stti Jc Jt 

/ / "^'^^"'^H + s)-H^w{T,x)dT. 



We can obtain operator convergence of the operator- valued integral in parenthe- 
ses if we deflect the contour C. Choose instead the portion C" of the hyperbola 
(See Figure [1]) 



(7) 



(where A is the smallest magnitude eigenvalue of —H) which lies in the left half- 
plane as our new contour. Then, since —H : C°'" C°^" is sectorial about 
(A — a)/2, Theorem 1.3.4 in [10] implies that the integral 

1 



2Tii 



s(r-t) 



{H + s)-'ds 



defines an operator- valued semigroup e ^^"^ *\ so the formula for F is given by 



(8) 



T{w){t,x) 



,-Hir-t] 



wIt, x)dT. 



It remains to show that the image of F is in fact Ya, as it is easy to see that its 
image is in L°° . That the image is as advertised is not immediately obvious because 
the contour deflection C ^ C changes the domain of the Laplace transform. In 
particular, the derivation given above is no longer valid with the new contour. 

Therefore, we must estimate ||f||za (recall that A is the smallest magnitude 
eigenvalue oi —H) 
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where < Ki, K2 < 00 are independent of t and w. We have made use of the usual 
estimate of the norm of {H + s)~^ : C°'" — > C*^'" when s is in the resolvent set of 
—H. In particular, note that the choice of C" being to the left of —a is crucial to the 
convergence of the integrals. Thus the image of F lies in Za- The backward-time 
decay of ^ is immediate from the Harnack inequality, so in fact the image of T lies 

inr,. 

Theorem 5. The linear map L : Ya{{— 00, 0]) — > oo,0]) is surjective and 

has a finite dimensional kernel. Therefore the set N~^{0) is a finite dimensional 
manifold, which is the unstable manifold of the equilibrium f . The dimension of 
N^^{0) is precisely the dimension of the positive eigenspace of H . 

Proof. The only thing which remains to be shown is that the domain Ya splits into 
a pair of closed complementary subspaces: the kernel of L and its complement. 
That its complement is closed follows immediately from a standard application of 
the Hahn-Banach theorem. (Extend idkcr l to all of Ya-) O 

Combining the fact that an equilibrium solution can have an empty unstable 
manifold (a numerical computation of the dimension of the eigenspaces of L can be 
found in [17] ) and is yet unstable, we have proven the following result. 

Theorem 6. All equilbrium solutions to ([l} are degenerate critical points in the 
sense of Morse. 

5. Linearization about heteroclinic orbits 

We can extend the technique of the previous section to the linearization about a 
heteroclinic orbit. The resulting generalization of Theorem [5] is that the connecting 
manifolds of H]) are all finite dimensional. 

Suppose that u is a heteroclinic orbit of ([1]) . Let /_ , /+ be the equilibrium 
solutions of ID) to which u converges as t —> —00 and t — > +00 respectively. 

Suppose that Aq : M ^ (0, 00) is the smallest positive eigenvalue of H{t). It is 
easy to see that Aq is piecewise C^, for instance, see Proposition 1.7.2 in [12 . The 
fact that the the spectrum of H lies entirely to the left of max{2||/+||oo, 2||/_||oo} 
ensures that Aq is a bounded function. We will define a pair of bounded, piecewise 
functions Ai and A2 which will aid us in defining a two more pairs of function 
spaces. Let Ai : M ^ (0, 00) be a bounded, piecewise function with bounded 
derivative which has the following properties: 
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Figure 2. Definition of Ai and A2 



Xi{t) is never an eigenvalue of H{t), 

Ai(t) 



limf^_ 
since u 



Ao(t) 



< 1, and 



f± uniformly, for a sufficiently large i? > 0, Ai can be chosen so 
that there are no jumps on its restriction to M — [—R, R]. 
Defining A2 is a somewhat more delicate problem. We would like to exclude 
the solutions which lie in the unstable manifold of /+, since they cannot lie in the 
space of heteroclines from /_ /+. We do this by separating the eigenvalues 
corresponding to the intersection of the unstable manifolds of /_ and /_(_ from 
those which lie in the stable manifold of f^. However, there is an obstruction to 
this technique. In particular, the eigenvalues of H{t) = — 2u{t) vary with time, 
and can bifurcate. To avoid this issue, we need some kind of regularity for the 
eigenvalues to prevent them from bifurcating. We follow Floer j7j in the following 
way: 

Conjecture 7. There is a generic subset (a Baire subset) of choices for in ([T|) so 
that if w is a heteroclinic orbit, all of the eigenvalues of H{t) are simple. 

Numerical evidence, as exhibited in [T7] and Section [5] suggests that the above 
Conjecture is true. When we assume that all of the eigenvalues of H{t) are simple, 
and therefore do not undergo any bifurcations other than passing through zero, we 
shall say u is a heterocline contained in Ureg- 

Let A2 be in (M) such that 

• A2 = Ai on [R, 00), and 

• X2{t) is not an eigenvalue of H{t) for any t. 
We can do this when u G Ureg- See Figure [21 

Definition 8. Define the Banach algebra Y\.{X) (for i = 1,2) to be the set of u 

in (X, C°'"(M)) such that the norm 



< 00, 
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where X is an interval containing zero. Likewise, we can define the spaces {X) C 
C*'(X, C'''"(M)) in a similar way. That these are Banach spaces follows from the 
boundedness of the A^. It is also elementary to see that these are Banach algebras. 

We then consider Ni , Li as Y\. (R) — ^ Zx^ (R) , where Li is the linearization of 
Ni about u for i = 1,2. (Again, since Y\. and Zx^ are Banach algebras, Li is the 
Frechet derivative of Ni.) For a i G {1,2}, consider the restriction of L,; to a 
map Fa, ((— cxo, 0]) Za, ((-oo, 0]). We rewrite 

(9) ^^^-U-a^ + 2/-j+(2/--2.). 

Likewise, we can define Lf : Yxj[0, oo)) Za;([0, oo)). 

We define the positive eigenspaces V+ for the equilibria as well 
(10) 

V+{f±) = span |t) e C°'"(R)| there is a A > with - 2/±^ v = Xv 

Note in particular that dim V'^{f±) < oo. 

Lemma 9. If u E Ureg is a heterocline that converges to f± as t ±oo, then the 
operator Li has a finite dimensional kernel for i G {1,2}, and in particular 

lim dim {u{t)) — lim dim {u{t)) < dim ker Li < dim ker L^ < oo. 

t — > — oo t — *+oo 

(The condition u G Ureg is only necessary for the i — 2 case.) 

Proof. Notice that the first term of ^ has finite dimensional kernel by Lemma |3] 
and closed image by Theorem[S] The second term of ([S]) is a compact operator since 
M ^ /_ uniformly. Thus L^ has a finite dimensional kernel. Let span{wm}^^j^ = 
ker L^ and consider the set of Cauchy problems 

f # = #4 - 2u/i for i > 
(11) <^ 9* ^ 

[h{0,x) = v„i{0,x). 

Standard parabolic theory gives uniqueness of solutions to (fTTj) , and that a solution 
h lies in the kernel of Lf, the restriction oi Li to [0, oo) xR. Therefore dim ker Li < 
dim ker L~ < oo. 

For the other inequality, modify u outside of [—R, i?] x R to get a u so that the 
linearization Li of N about u satisfies 

• ker Li is isomorphic to ker Li as vector spaces, 

• w|(-oo,-fl)xR = ./-, and 

• w|(fi,oo)xR = /+• 

We can do this for a sufficiently large R, since u tends uniformly to equilibria. Then 
the flow of 

dh _ d^h 
dt dx^ 

defines an injective linear map from the timeslice at —R to the timeslice at R. (That 
is, it gives an injective map from C'^'"(R) to itself - injectivity being an expression 
of the uniqueness of solutions.) Each clement v of the kernel of Li evidently must 
have v{—R) G V'^{f-) and v{R) ^ V~^{f+). Therefore, the injectivity ensures that 
the intersection of the image under the flow of V~^{f-) with the complement of 
V^{f+) has at least dimension dim V^{f-) — dim F+(/+). □ 
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Remark 10. Multiplication by u, C^R^, C°'"(M)) ^ C°{R'^) is not a compact 
operator, in particular note that dim ker = oo. 

Theorem 11. Let u be a heterocline of ^ which connects equilibria f±. There 
exists a union [J Mu of finite dimensional submanifolds ofC^{R, C'''"(R)) which 

• contains u and 

• consists of heteroclines connecting /_ to /_|_ . 

Ifu e Ureg, then has dimensionYivcit^-oo dim {u{t))—\mit^oo dim V'^{u{ty), 
and this is maximal among such submanifolds M„. 

Proof. Observe that Li is surjective, since it is easy to show that the formula 

Ti{w){t)^ e-i'»''"^^^'^^w{T,x)dT 
is a well defined right inverse of Li. This involves showing that 

27ri Jc{t} 

converges, where we note that the contour changes with time. As it happens, the 
computation in [10] goes through with the only change that at t = 0, we deflect the 
contour to the right, rather than the left (as in Figure [1]). Since Lemma [9] shows 
that Li has finite dimensional kernel, then it follows that M„ = N^^{0) is a union 
of finite dimensional manifolds, with a finite maximal dimension. It is obvious that 
Mu consists entirely of heteroclinic orbits and contains u. 

It remains to show that the dimension of Af„ is as advertised and maximal. 
Observe that L2 is a compact perturbation of an operator L'2 : i^) (R) 
which is time-translation invariant. This follows from the precise choice of A2 being 
continous and not intersecting the eigenvalues of H. L2 and L'2 are both surjective 
by exactly the same reasoning as for Li. L^ is injective by using separation of 
variables as in Lemma |3] (noting that all nontrivial solutions blow up in the Y\.^ 
norm). Therefore the Fredholm index of L'2, hence L2 is zero. However, this implies 
that L2 is injective. 

Since L2 is bijective, any solution to L2U = which decays faster than ^'^'^^^'^^ 
as t —00 ends up growing faster than e-/^^^ as t ^ +00, and in particular 
does not tend to zero. As a result, such a solution cannot be in ker Li. This implies 
that dim ker Li < limt^_oo dim V~^{u{t)) — linit^oo dim V^{u{t)), which with the 
estimate in Lemma [9] completes the proof. □ 

Remark 12. Even if w ^ Ureg (when there exist nonsimple eigenvalues of H{t)), the 
function Ai can still be constructed. As a result, we always get that the connecting 
manifold Mu is finite-dimensional. 

Corollary 13. The space of heteroclinic orbits has the structure of a cell com- 
plex with finite dimensional cells. This cell complex structure is evidently finite 
dimensional if there exist only finitely many equilibria for ([T]) . 

6. An extended example 
Consider the following spectial case of H]) 



(12) 



du d^u , , , . _ 
-dt-dx^-'' 
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where the choice of in ([T]) has been fixed. The bifurcation diagram for the 
equihbria of (|12p can be found in Figure[31 The bifurcation diagram is parametrized 
by three variables: c, /(O), /'(O). (Since the equihbrium equation is a second-order 
ODE, it suffices to specify each solution by its value and first derivative at 0.) 
Based on the Theorem [Sj the number of positive eigenvalues shown in Figure [3] 
corresponds exactly to the dimension of the unstable manifold of each equilibrium. 
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Figure 4. Behavior of solutions near the frontier of the stable 
manifold of /o (horizontal axis is x) 

6.1. Frontier of the stable manifold. According to Figure [3l when c = —1.2, 
there is only one equilibrium, /q. It has empty unstable manifold, though of course 
it is asymptotically unstable (as is shown in [l9|). On the other hand, /o has an 
infinite dimensional stable manifold, which is not all of C°'"(K), as a consequence of 
the asymptotic instability. As a result, its stable manifold has a frontier in C°'"(R) 
(which may not be a boundary in the sense of a manifold with boundary) . We are 
interested in the qualitative behavior of solutions near and along this frontier. We 
know by Lemma 6 of [18] that if they tend to /o uniformly on compact subsets, 
then they do so uniformly. It is enlightening to use a numerical procedure to this 
end. We start solutions at the following family of initial conditions 

(13) ua(x)^ h{x)+Ae-^"'^''. 

Using the Fujita technique (exactly as shown in [19]), we can show that for suffi- 
ciently negative the solution started at ua will not be eternal. As a result, the 
family of initial conditions ua intersects the frontier of the stable manifold of /q. 
An approximation to the value of A which corresponds to the frontier can be easily 
found using a binary search. Some typical such solutions are shown in Figure U] 
and the approximate value of A corresponding to the frontier is A « —2.15 

The qualitative behavior shown in Figure |4] indicates that there is some kind of 
traveling disturbance in the frontier solutions, which seems like a traveling wave. 
However, such a solution also appears to tend uniformly on compact subsets to /o, 
so in fact it converges uniformly. (The uniform convergence is not obvious from the 
figure, due to the numerical solution being truncated at a finite time.) The leading 
edge of this disturbance collapses to ~oo in finite time for solutions just outside 
the stable manifold of /q. 

6.2. Flow near equilibria Virith tv^ro-dimensional unstable manifolds. Also 
of interest is the structure of the flow in the unstable manifold of the "fork arms" 
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Figure 5. Flow in the unstable manifold of a "fork arm." c = 
0.0600 (left); c = 0.0501 (right) 

which occur a,t c — 0.0740, as they approach the pitchfork bifurcation at c = 0.0501. 
Figure [5] shows a schematic of the flow based on numerical evidence. Of particular 
interest is the behavior near the boundary marked A. Solutions to the right of the 
boundary are not eternal solutions - they fail to exist for all t. Solutions to the left 
of A are heteroclinic orbits connecting the equilibrium with an unstable manifold 
of dimension 2 to the equilibrium with an unstable manifold of dimension zero. A 
typical such solution is shown in Figure [H 

To examine solutions near the boundary A, we center our attention on the case 
c = 0, which has two equilibria, one of which (call it /i) has a 2-dimensional 
unstable manifold. (This corresponds to the right pane of FigureO) If we linearize 
about /i, the operator H ^ - 2fi : C°^"(M) ^ C°'"(M) has a pair of simple 
eigenvalues, as is easily seen in the right pane of Figure [H] at t — 0. One of these 
eigenvalues is smaller, to which is associated the eigenfunction ei in Figure [71 The 
eigenfunction 62 is associated to the larger eigenvalue. In Figure El ei corresponds 
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Figure 7. Eigenfunctions describing unstable directions at /i 

to the horizontal direction, and 62 corresponds to the vertical direction. From the 
proof of Lemma [4l it is clear that {61,62} spans the tangent space of the unstable 
manifold at /i. Therefore, we specify initial conditions UA.eix) for a numerical 
solver using 

(14) UA,e{x) = fi{x) + A{ei{x) cosO + e2{x)sm.e) . 

(Taking A small allows us to approximate solutions which tend to /i in back- 
wards time.) Since the perturbations along 61,62 are quite small, and indeed the 
eigenvalue associated to ei is much smaller than that associated to 62, examining 
the numerical results of evolving UA,e is quite difficult. The behavior along the 
boundary occurs at a much smaller scale than /i, yet is crucial in determining the 
long-time behavior of the solution. To remedy this, the boundary behavior is bet- 
ter emphasized by plotting UA,e{t, x) — fi{x) instead. Figure [5] shows the results of 
evolving initial conditions p4| for A — OA and various values of 6. 

Solutions in Figure [8] show a similar kind of behavior as in the case of the 
frontier of /q. There is a traveling front, which moves very slowly in the negative 
cc-direction. However, the behavior is quite a bit more delicate. The determining 
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0=1.11494 



0=1.11496 



0=1.11497 




Figure 8. Difference between equilibrium /i and the numerical 
solution started at UA,e, where black indicates a value of -0.2, and 
white indicates 0.2. The horizontal axis represents t, and the ver- 
tical axis represents x. A = 0.1 in all figures. Starting from the 
upper left, = 1.11494, 1.11496, 1.11497, 1.11498, 1.11499, 1.115. 



factor in locating the frontier of /o is the perturbation in a direction roughly like 
62, which has a large eigenvalue. On the other hand, for /i. Figure [5] indicates that 
such a direction is not parallel to the boundary of the connecting manifold. (The 
boundary direction is some linear combination of ei and 62, with a numerical value 
for the angle 9 being roughly 1.114975 radians.) The eigenvalue associated to ei is 
roughly ten times smaller, and therefore perturbations in that direction are much 
more sensitive. Additionally, the action of the flow is therefore primarily in the 
direction of ei, which tends to mask effects in other directions. For this reason, 
it was visually necessary to postprocess the numerical solutions by subtracting /i 
from them. Otherwise the presence of the traveling front was unclear. 
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7. Conclusions 

We have shown that the tangent space at an equihbrium sphts into a finite di- 
mensional unstable subspace, and infinite dimensional center and stable subspaces. 
However, it is quite clear by [19j that the center subspace is nonempty and large. 
Indeed, considering the work of [20], the center and stable subspaces are not closed 
complements of each other. Additionally, we have given conditions for the space of 
heteroclinic orbits to have a finite dimensional cell complex structure. 
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